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Relative Motion Using Analytical Differential Gravity 


1.0 INTRODUCTION 

This paper presents a new approach to the computation of the motion of one satellite 
relative to another. The trajectory of the reference satellite is computed accurately 
subject to geopotential perturbations. This precise trajectory is used as a reference 
in computing the position of a nearby body, or bodies. 

The problem that arises in this approach is differencing nearly equal terms in the 
geopotential model, especially as the separation of the reference and nearby bodies 
approaches zero. By developing closed form expressions for differences in higher 
order and degree geopotential terms, the numerical problem inherent in the 
differencing approach is eliminated. 
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2.0 ANALYSIS 


The equations of motion for a satellite moving under the influence of gravity are 
written 


av 



where r T = ( ri r 2 r 3 ) 

where V is the potential function 


( 1 ) 


V = - 



E -B- 

m =0 r 



(Cnm cos m A + Snm sin m A) 


( 2 ) 


and where P n m(£) is the Legendre polynomial of degree n and order m, £ is r 3 /r, C n m 
and Snm are the geopotential model coefficients, tan A = r 2 /r v and a e is earth 
equatorial radius. 


Using the recursive formulation given in [2], and considering only terms through 

• • 

n = 2, m = 2, the equation for r may be written: 


• • 

r = - 


a e 
r , 


{(r 2 + EH 2 )-M k 2 )} 


(3) 


where r 2t J 2( K 2 an d H 2 may be shown to be 

r 2 = — G — 020 

r r 2 


where 


G = 


15C 

0 
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(12S 21 + 30S 22 ) 
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where 

f = (6C 22 6S 22 3C 21 ) 

K 2 = k T 

2 - r 

where 

k T = (6S 22 6C 22 3C 21 ) 


and 


r 



where 

h T = (3C 21 3S 21 3C 20 ) 


Using these, we may write equation (3) as, 



by noting that r 3 may be written 
r a =r T a 

where 

a = 
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and then defining 


0 — G + a h 


0 = 1 - 


k T 



we may write equation (4) as 


** r 3 ~ 2 r 2 r 

r = - p O -=- H — P '-'20 a e — - p a e ( r T 0 r ) -=- 


(5) 


In fact, it can be shown that in general for the geopotential 


• • 

r ss Mr 

where M is a matrix whose elements may depend on r. 


Assume that we have two satellites, the first with position vector r and the second 
with position vector p. Both must satisfy equation (5). Subtracting the two 

differential equations and defining the difference between the two solution vectors 
to be 

8 = p— r 


we may write 


8 = — p <J> 



C 2 
+ 3 p a e 

2 




2 
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(p T 0p)JL 

“ “ P 


~ ~ -7 


( 6 ) 
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Now collect coefficients of — * — » to get 

r 3 r 5 r 7 


-P 


20 a e 


8= — 4»(8+( 1)_p ) + 3p — — l 8 + (— -l)p 


2 r 5 


-pa e J 


(p T 0p)( — -1) +2 (8 T 0p)-(8 T 08) 


P 7 


p -l-(r T 0r) 8 


1 


(7) 


Note that the factors 


— -1 , — -1 , and — -1 
P 3 P 5 P 7 

should each approach zero as 8 approaches zero. Numerically this presents a 
problem since r and p are large and nearly equal. 

These can all be computed using Potter’s 111 approach by noting that 


r 2 = (p-8) . ( p— 8) = p 2 - 2 p • 8 +8 2 


and 


r"= ( p 2 — 2 p • 8 +8 2 ) 


n/2 


and that 


p" = (p 2 ) 


n/2 


p 2 - 2 p • 8 + 8 2 


n/2 


-1 = 


- 1 


( 8 ) 
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and we can write 


£l ! _ ( 1 + q) n - 1 

p n (1+q) n/2 +1 


(9) 


where 


8 2 - 2 p . 8 

q a = 


Note that q — 0 as 8 — ► 0 


We can now write these factors as 


r 3 


P 


3 


- 1 


f 

(l + q) 3/2 + 1 


where 


f = 3q + 3q 2 + q 3 


and 




- 1 = 


f + (f+1)(2q+ q 2 ) 
(1+q) 5/2 + 1 

q + (q+1)(2f+f 2 ) 
(1+q) 7/2 + 1 


( 10 ) 


(ID 


( 12 ) 


Note that from the definition of q and f, these factors approach zero directly as 8 
approaches zero. 
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3.0 DISCUSSION AND RECOMMENDATIONS 


Equations (10), (11), (12) when substituted into equation (7) yield the companion set 
of differential equations for the second satellite relative to the first. The terms in the 
resulting equation all go to zero directly as 8 approaches zero and do not contain 
differences of large nearly equal terms. These equations would be quite useful for 
both space station and tethered satellite analysis. 

The technique presented here extends, albeit with effort, to higher order and degree 
terms in the geopotential. A recursive approach to the computation of the coefficient 
matrices would be a welcome development. 
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